Detection of stacking faults breaking the [110]/ [110] symmetry in ferromagnetic 
semiconductors (Ga,Mn)As and (Ga,Mn)(As,P) 
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We report high resolution x-ray diffraction measurements of (Ga,Mn)As and (Ga,Mn)(As,P) epi- 
layers. We observe a structural anisotropy in the form of stacking faults which are present in the 
(111) and (111) planes and absent in the (111) and (111) planes. The stacking faults produce no 
macroscopic strain. They occupy 1(T 2 - KT 1 per cent of the epilayer volume. Full-potential den- 
sity functional calculations evidence an attraction of Mnc a impurities to the stacking faults. We 
argue that the enhanced Mn density along the common [110] direction of the stacking fault planes 
produces sufficiently strong [110]/ [110] symmetry breaking mechanism to account for the in-plane 
uniaxial magnetocrystalline anisotropy of these ferromagnetic semiconductors. 

PACS numbers: 61.72.Dd, 75.50.Pp, 78.55.Cr 



The rich phenomenology of magnetocrystalline 
anisotropics in ferromagnetic (III,Mn)V semiconduc- 
tors has been the prerequisite of numerous studies 
of magnetic, magneto-transport, and magneto-optical 
phenomena and of prototype semiconductor spintronic 
devices pQ. The in-plane biaxial anisotropy as well 
as the out-of-plane uniaxial anisotropy terms are well 
understood based on the periodic crystal structure char- 
acteristics of the epilayers [2 j. The former term is due 
to the cubic symmetry of the host III-V semiconductor 
and the latter term due to the lattice-matching strain 
induced by the difference between lattice parameters of 
the free-standing (III,Mn)V crystal and the substrate. 

The most extensively studied material is (Ga,Mn)As 
grown on GaAs, in which the compressive strain in 
the ferromagnetic (Ga,Mn)As epilayer makes the out-of- 
plane orientation of magnetization energetically unfavor- 
able. It has been recognized from the early studies of 
these in-plane ferromagnets that the biaxial anisotropy 
term is complemented by an additional uniaxial term 
breaking the symmetry between [110] and [110] crystal 
directions. While the presence of this term and its com- 
petition with the biaxial anisotropy term have played a 
key role in the research of ferromagnetic semiconductors, 
including studies of electrical or optical manipulation of 
the magnetic state PQ, the microscopic origin of the in- 
plane uniaxial anisotropy has remained elusive. 

We have considered the following guidelines when 
searching for the uniaxial in-plane symmetry breaking 
mechanism: (i) The corresponding uniaxial magnetocrys- 
talline anisotropy is not a surface or interface effect but 
a bulk phenomenon, (ii) No strain component has been 
detected in the epilayers that would break the symme- 
try between the two in-plane diagonals, (iii) No system- 
atic dependence has been identified in the in-plane uni- 
axial magnetic anisotropy on the growth induced lattice- 
matching strain, as seen e.g. from the comparison of com- 



pressively strained (Ga,Mn)As/GaAs epilayers and ten- 
sile strained (Ga,Mn)(As,P)/GaAs epilayers. (iv) From 
the effective modeling of the in-plane uniaxial anisotropy 
it has been concluded that the symmetry breaking mech- 
anism is related to the high MnQ a doping [2 . 

In this paper we report experimental observation and 
theoretical investigation of stacking faults in (Ga,Mn)As 
and (Ga,Mn)(As,P) which break the in-plane [110]/ [110] 
symmetry and whose characteristics are consistent with 
the above guidelines. The experimentally estimated 
density of the stacking faults and the theoretically in- 
ferred attraction of MnQ a to these lattice defects yields 
a strength and sense of the symmetry breaking mecha- 
nism which we compare to the broken crystal symmetry 
due to an in-plane uniaxial strain. The latter symme- 
try breaking mechanism has been commonly used as an 
effective parameter to theoretically model the uniaxial 
magnetocrystalline energy in unpatterned epilayers or as 
a real tool to control the anisotropy in microstructured 
films or in epilayers attached to piezostressors [2 . 

Measured (Ga,Mn)As and (Ga,Mn)(As,P) samples 
were grown by low-temperature (200-230° C) molecu- 
lar beam epitaxy (LT-MBE) on a GaAs substrate and 
buffer layer. For more details on the sample growth see 
Refs. [3[ H]. X-ray experiments were carried out on the 
diffraction beamline at the ELETTRA synchrotron fa- 
cility in Trieste. Samples were mounted on a two-axis 
tilt platform and aligned in the way that the normal to 
the sample surface coincides with the rotation axis of the 
diffractometer; the rotation angle is (j). In this geometry, 
the glancing angle between the beam and the sample 
surface remained constant during the 0-scan. The energy 
of the incident beam was set to 10.3 keV, i.e. just below 
the K absorption edge of gallium, in order to minimize 
the absorption in the sample and to avoid fluorescence. 
The beam size was set to 500 /am (horizontal) and 20 /im 
(vertical) using two pairs of slits in front of the sample. 
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A small vertical aperture guarantees the elimination of 
the undesired scattering from sample borders. The area 
detector Pilatus 2M (1475 x 1679 pixels of the size of 
172 x 172 jim 2 in 3 x 8 modules, dynamic range of 20 
bits) made it possible to collect the weak diffuse scatter- 
ing patterns using the exposure rate of 5 s per frame. One 
scan was composed of 720 frames collected with the an- 
gular step Acj) = 0.5°. The frames were put together and 
the intensity maps in the large volume of the reciprocal 
space were constructed. 

In Figs. l(a),(b) we show measured data for 500 nm 
thick 5% Mn-doped (Ga,Mn)As epilayer. The large sam- 
ple thickness allows us to set the glancing angle 6 to 0.3°, 
which is slightly above the critical angle 6 C = 0.25°. Un- 
der these conditions a large portion of the (Ga,Mn)As 
epilayer volume is illuminated while only a small por- 
tion of the beam penetrates into the substrate and the 
scattering from the substrate is negligible. Two cross- 
sections through the reciprocal space perpendicular to 
each other are shown in Figs. l(a),(b). The cross-section 
in panel 1(a) contains only diffraction spots along with 
truncation rods perpendicular to the sample surface. On 
the other hand, the cross-section shown in panel 1(b) 
contains additional diffuse streaks in directions [111] and 
[111], indicating the presence of stacking faults in planes 
(111) and (111). 

In Fig. 2(a) we present a detailed image of the area 
near the -111 diffraction and compare in Fig. 2(b) with 
the same area measured at = 0.2°, i.e., slightly be- 
low the critical angle. In the latter image the thick- 
ness of the illuminated (Ga,Mn)As surface layer is only 
~ 10 nm. Consistently, the truncation rod perpendicular 
to the sample surface is more extended in Fig. 2(b) than 
2(a). The comparison of the two images demonstrates 
that the stacking faults occupy a significant part of the 
epilayer volume since the intensity of the diffuse [111] and 
[111] streaks is significantly larger in Fig. 2(a) than 2(b). 

We have performed additional measurements at = 
0.2° on (Ga,Mn)As and (Ga,Mn)(As,P) epilayers of 
width ranging from 35 to 100 nm, as-grown and annealed, 
and with Mn-doping up to 10% and P-doping of 9%. In 
all cases we obtained results similar to Fig. 2(b) with the 
diffuse streaks only in directions [111] and [111]. No dif- 
fuse streaks in any direction were observed in a reference, 
LT-MBE grown undoped 500 nm thick GaAs epilayers. 
From these data we can conclude that the stacking faults 
span across the epilayer volume, are related to the Mn- 
doping, are present in the same planes in samples with 
smaller (annealed GaMnAs) and larger (as-grown GaM- 
nAs) magnitude of the growth lattice-matching strain 
and in epilayers with both compressive (GaMnAs) and 
tensile (GaMnAsP) strain. We also point out that the 
stacking faults are present in thick epilayers as well as in 
thin, annealed high magnetic quality films. 

From the modeling of the measured x-ray scattering 
patters (Figs. 1(c) and 2(c)), we identified the micro- 
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FIG. 1: (a, b) Cuts through the reciprocal space of the mea- 
sured diffraction intensity maps. The vertical truncation rods 
perpendicular to the sample surface are due to finite thickness 
of the illuminated film. The streaks in directions [111] and 
[11-1] represent the diffraction on stacking faults in planes 
(111) and (11-1), resp. (c)The calculated x-ray diffraction 
considering the distribution of the stacking faults described 
in the text for the cut through the reciprocal space as in (b). 



scopic nature of the stacking faults. In the unperturbed 
GaAs zinc-blende crystal structure, the plane stacking 
along the [111] direction (or any other body diagonal) 
is characterized by a repeating a-b-c sequence of Ga-As 
planes, as shown in Fig. 3(a). A stacking fault along the 
[111] diagonal shown in Fig. 3(b) corresponds to one of 
the planes missing in the stacking. The fault can be also 
viewed as replacing the zinc-blende stacking by two in- 
terpenetrating wurtzite stackings highlighted by the two 
red (vertical) bars in Fig. 3(b). A stacking fault due 
to an extra plane, or more separated two wurtzite se- 
quences, is depicted in Fig. 3(c). A removed and im- 
mediately inserted plane creates a fault comprising two 
neighboring wurtzite stackings, as shown in Fig. 3(d). 
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FIG. 2: Detail of the diffuse scattering in the vicinity of the 
reflection -111 measured at (a) = 0.3° and (b) = 0.2°. (c) 
Measured intensity profile along the [111] direction passing 
through the reflection -111 and corresponding numerical sim- 
ulation. Here q x is the x co-ordinate of the diffraction vector 
and a* is the reciprocal lattice unit cell vector. 



Our experimental data could not be fitted by the stack- 
ing faults in 3(b) and 3(c). As seen from the compari- 
son with Fig. 3(a), these two faults yield a zinc-blende 
stacking after the fault which is shifted by one plane with 
respect to the unperturbed lattice. It disrupts the match- 
ing of the epilayer with these stacking faults to the sub- 
strate on a macroscopic scale and, therefore, the faults 
are unlikely to form. Their absence is also consistent 
with the absence of an in-plane uniaxial strain which we 
confirmed experimentally in our epilayers with accuracy 
~ 10 -5 . (We compared average lattice parameters in the 
[110] and [110] directions by using relative position of 17 
Bragg reflections of the Mn-doped epilayer with respect 
to these of the GaAs substrate.) The stacking fault in 
Fig. 3(d), on the other hand, produces no macroscopic 
strain and indeed our fittings show that this fault repre- 
sents ~ 90% of the stacking faults present in the mea- 
sured epilayers. The remaining faults, producing also no 
macroscopic strain outside the fault, have one or more 
zinc-blende a-b-c sequences between the removed and in- 
serted plane and their probability with increasing num- 
ber of inner a-b-c sequences was fitted by a T-function 
distribution. 

As pointed out in the previous paragraph, there is no 
experimentally observable strain in the (Ga,Mn)As epi- 
layers, within the ~ 10 -5 error bar which would break 
the [110]/ [110] crystal symmetry. The uniaxial in-plane 
strain has been used as a hypothetical symmetry break- 
ing mechanism to account in microscopic theory calcula- 
tions for the experimentally observed uniaxial in-plane 
magneto-crystalline anisotropy [2j |6]. Typical magni- 
tudes of the strain required by the modeling to reproduce 
the measured magnitude of the magnetic anisotropy are 
~ 10 -4 , i.e., an order of magnitude larger than the ex- 
perimental error of our measurements. Based on this 
we can exclude strain as the real microstructural origin 
of the [110]/ [110] magnetic anisotropy. Results of our 
ab initio calculations, described in the following para- 
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FIG. 3: (a) Zinc-blende stacking in an ideal GaAs lattice. 
(b)-(d) Stacking faults with one removed, added, or removed 
and immediately inserted planes, resp. (e) Calculations of 
the total energy of the supercell with the respective stacking 
fault as a function of the position of the substitutional Mnc a . 
Total energies are plotted with respect to the energy of the 
supercell with Mno a at the 0-th layer in the supercell. 



graphs, suggest that the observed stacking faults pro- 
vide the searched for microstructural symmetry breaking 
mechanism. Our supercell calculations were performed 
using the full-potential linearized- augmented-plane- wave 
method with the improved local-density approximation 
by the generalized-gradient approximation (LAPW-GGA 
WIEN package [5]). The calculated formation energies 
of the stacking faults shown in Figs. 3(b)- (d) for an un- 
doped GaAs are in the range of 13-15 mRy. In Fig. 3(e) 
we plot the dependence of the formation energy of the 
stacking faults on the position of the substitutional MnQ a 
impurity along the [111] (or [111]) crystal direction in the 
supercell. By comparing the results for the three consid- 
ered stacking faults (Figs. 3(b)-(d)) we can conclude that 
MnQa is attracted to the wurtzite stacking perturbation 
in the lattice and the energy profile can be understood 
as a convolution of attractive potentials due to the indi- 
vidual wurtzite perturbations. Remarkably, in our super- 
cells corresponding to 25% of Mnc a impurities, the depth 
of the attractive potential energy well is as large as 20% 
of the formation energy of the stacking fault in GaAs. 
We conclude that Mn impurities have a strong tendency 
to decorate the stacking faults and that by this process 
the Mn impurities significantly enhance the probability 
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FIG. 4: (a) Schematic image of uniformly randomly dis- 
tributed Mn (red (dark) dots) in the epilayer plane in an 
unstrained lattice with cubic symmetry, (b) Schematics of 
the planes in which the stacking faults are observed by x-ray 
diffraction, (c) Schematic image of anisotropic Mn distribu- 
tion in the epilayer plane in an unstrained lattice which results 
from the tendency of Mn to decorate the stacking faults, (d) 
Schematic image of uniformly distributed Mn on the lattice 
sites in the presence of a uniaxial in-plane symmetry breaking 
strain, (c) and (d) are analogous from the perspective of the 
anisotropic mean separation of Mn in the epilayer plane. 



of the formation of the stacking faults in the epilayer. 

In schematic diagrams in Fig. 4 we illustrate the in- 
ferred connection between the Mnc a decorated stacking 
faults and the uniaxial magnetocrystalline anisotropy. 
Fig. 4(a) shows schematically the biaxial symmetry in 
case of unperturbed lattice and uniform distribution of 
MnQa impurities. Since the stacking faults are observed 
only in the (111) and (111) planes (Fig. 4(b)) and Mn Ga 
is attracted to the stacking faults, we can expect higher 
mean density of the MnQa impurities along the com- 
mon [110] direction of the two planes. This leads to an 
anisotropic distribution of MnQa impurities in the epi- 
layer plane with a smaller mean distance between impu- 
rities in the [110] direction than [110] direction. 

In Fig. 4(d) we illustrate that the non-uniform Mnc a 
distribution in an unstrained square lattice is from the 
perspective of the mean MnQa separation analogous to 
the case of a uniform Mn distribution and a strained 
lattice. The latter case has been considered previously 
in the effective modelings of the corresponding uniaxial 
magnetic anisotropy. Since the density of the stacking 
faults we observe is ~ 0.01 — 0.1%, the corresponding 
symmetry breaking mechanism is strong enough to ex- 
plain the uniaxial magnetic anisotropy. (Compare the 
density of the stacking faults with the ~ 10 -4 strain used 
in the effective modeling of the magnetic anisotropy.) We 
note that detection of the stacking faults and of the in- 
creased MnGa density on these defects is likely unachiev- 
able by direct atomic resolution scanning techniques. 
This underlines the merit of our approach, combining 
the macro-scale x-ray measurements with ab initio total 



energy calculations, for identifying the subtle symmetry 
breaking mechanism. 

The sign of the magnetic anisotropy in the Mn-doped 
magnetic epilayers is not uniquely determined by the sign 
of the underlying broken structural symmetry; it can vary 
for the given sense of the broken symmetry as a function 
of, e.g., the hole density in the magnetic film. Typi- 
cally, however, the [110] direction is magnetically eas- 
ier than [110] direction. The common [110] direction of 
the stacking fault planes then coincides with the typical 
easy magnetic direction. If again considering the anal- 
ogy with the case of uniform Mnc a distribution and in- 
plane strained lattice, this observation provides further 
support for the stacking-fault origin of the uniaxial in- 
plane anisotropy. From the modelings considering the 
in-plane strain, as well as from experiments on intention- 
ally in-plane strained (Ga,Mn)As epilayers by lithograph- 
ical pattering or piezostressors [2] , it has been concluded 
that the magnetic easy direction is along the smaller 
lattice parameter, i.e., along the direction with smaller 
mean separation between Mnc a impurities. Our results 
show that in unstrained epilayers the stacking faults yield 
smaller mean separation between Mnc a impurities along 
the in-plane axis which also coincides with the typical 
uniaxial magnetic easy direction. We conclude that all 
characteristics of the experimentally observed stacking 
faults in the studied (III,Mn)V epilayers are consistent 
with the guidelines (i)-(iv) summarized in the introduc- 
tion for identifying the microstructural [110]/ [110] sym- 
metry breaking mechanism in these magnetic semicon- 
ductors. 
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